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The  propagation  of  mechanical  disturbances  through 
solid  or  fluid  media  is  a  nonlinear  phenomenon  whenever 
the  disturbance,  or  wave,  has  a  finite  amplitude.  The 
term  finite  amplitude  refers  to  any  disturbance  where  the 
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.splacemer.t  is  not  inf initesmally  small  and 


where  the  superpos ition  principle  no  longer  holds.  One 
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~v.c  waves  present  in  one  section  oi  volume  01  tne  neciur. 

v.  ill  affect  one  another,  i .  e  .  ,  there  will  be  an  interaction 

between  them.  Two  such  interactions  are  the  three -phonon 

interaction  and  harmonic  generation. 

In  the  three-phonon  interaction,  two  noncollinear 

finite  amplitude  waves,  which  we  shall  call  the  primary 

waves ,  intersect  in  such  directions  that  xhey  interact  to 

produce  a  third,  or  resultant  wave.  Hollins  and  Taylor^" 

have  observed  three-phonon  interactions  for  bulk  waves  in 

solids  which  were  chosen  large  enough  so  that  boundary 

2 

effects  could  be  neglected.  Tanski  has  extended  the 
concept  to  three-phonon  interaction  of  surface  waves  on 
l_crystals  with  one  defined  boundary.  A  three-phonon  _ 
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interaction  was  detected  :y  Ire wer^  or.  a r.  isotropic  plate  • 
with  two  defined  boundaries . 

The  analytical  sclutior.  tc  the  prc tier,  cf  the  non¬ 
linear  three-phonon  interaction  in  a  bulk  isotropic  medium 
has  been  obtained  by  Jones  and  Kobett  by  a  perturbative 

approach  and  Jreen's  function  techniques.  For  the  surface 

.  2 

wave  three-phonon  interaction  on  crystals,  Tarski  has 
suggested  several  possible  approaches  to  obtain  an  appro¬ 
ximate  solution.  The  first  goal  of  the  present  work  is 
to  find  the  solutions  for  the  three-phonon  interaction  in 
an  isotropic  plate  with  two  boundary  conditions. 

In  the  harmonic  generation  interaction,  specifically 
second -harmonic  generation,  an  acoustic  wave  generates 
higher  harmonics  while  propagating  in  a  solid  medium. 
Second-harmonic  generation  has  been  experimentally  observed 
for  bulk  waves  in  an  isotropic  medium  by  Viktorov^,  and 

Ldpen^1  has  studied  second-harmonic  generation  of  acoustic 

3 

surface  waves  on  crystalline  solids.  Erower  also  observed 
the  second-harmonic  generation  in  Lamb  modes  in  an  iso¬ 
tropic  plate. 

Different  solutions  to  the  problem  of  second-harmonic 
generation  in  a  bulk  isotropic  medium  have  been  obtained  by 
Kikata"  and  Viktorov-^.  The  second-harmonic  generation 
ijproblem  on  crystals  has  been  solved  approximately  by 


leper.',  who  used  a  standard  perturbative  technique.  The 
second  goal  of  this  thesis  is  to  find  the  solution  for 
second -harmonic  generation  in  an  isotropic  plate. 

In  this  study  a  theoretical  investigation  of  the 
acoustical  nonlinear  effects  for  an  isotropic  plate  will 
be  presented.  This  work  will  take  a  new  approach  for  the 
solution  to  the  nonlinear  equations  of  motion  for  an  iso¬ 
tropic  plate.  The  approach  involves  a  perturbative  appro¬ 
ximation  to  the  nonlinear  equations  of  metier,  and  a  sub¬ 
sequent  solution  by  the  Green's  function  technique.  The 
Ireer. '?  function  is  constructed  via  an  eigenfunction  ex- 
nans  ion.  Chatter  II  contains  a  review  of  the  three  -phono  r. 
interaction  and  second -harmonic  generation  m  a  bulk  medium 
and  on  a  crystal  surface.  Chapter  III  has  three  parts,  the 
first  of  which  is  a  review  of  Green’s  function  technique. 

In  the  second  part  the  orthogonal izat ion  relation  of 
eigenfunctions  for  a  plate  will  be  derived.  The  last 
..  ctior.  contains  the  derivations  leading  to  a  formal 
solution  for  the  nonlinear  equations  of  an  isotropic  plate. 
This  solution  for  the  special  cases  of  the  three-phonon 
interaction  and  second-harmonic  generation  will  be  inves¬ 
tigated  in  Chapter  IV  along  with  bounded  wave  considera¬ 
tions.  The  last  chapter  contains  the  conclusions  and 
for  fumure  worx. 
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CHAPTER  II 

THEORY 

A .  Acoustic  Nonlinear  Equations  of  motion 

The  equations  describing  the  motion  of  the  particle 
in  an  isotropic  solid  may  be  found  in  many  textDCOKS x’r . 
The  general  equation  of  motion  is 
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where  fQ  is  the  density  of  the  solid;  CiJc2ffl  and  CiJcpqrg 
are  the  second  and  third  order  elastic  constants,  respec¬ 
tively.  The  auantity  W  is  the  strain  tensor  and  may  be 
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Here  U  is  the  displacement  vector.  The  general  summation 
rule  of  repeated  indices  is  used  here  and  in  all  that 
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isotropic  solid,  where 
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■ere  K  is  the  compression  modulus,  u  is  the  snear  modulus, 
i.r.d  A,  3,  and  C  are  third  order  elastic  constants. 

In  practical  applications  a  perturbative  approach  to 
:r.is  ecuation  is  used  in  which  U  is  written  as 


U  =  U  +  U' 
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Here  UQ  is  the  solution  with  the  right  hand  side  of  Eq.II- 
3  equal  to  zero,  i.e.,  the  linear  homogeneous  equation. 

■  V 

Thus,  when  we  substitute  Eq.II-h  into  Eq.II-3»  UQ  dis¬ 
appears  from  the  left  hand  side  leaving  only  U'  on  the 
left.  On  the  right  hand  side  since  U’  is  small  compared 
to  UQ ,  it  is  neglected  from  terms  on  the  right.  So  the 
right  hand  side  of  Eq.II-3  contains  terms  in  UQ  only. 

I  Now,  the  nonlinear  partial  differential  equations  are 
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vector  F  is  acting  as  the  source  term  for  the 


wave.  ?  is  determined  by  substituting  U  in 


the 


side  of  £ 
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r .  Solving  the  Acousmic  Ecuatic ns  f:r  th>-  T:.ree  -iLr.or. 
Interac  x  lor. 

In  this  section  three-pi.cr.cn  interact ic;.s  in  isctrc 
bulk  media  and  on  crystalline  surfaces  will  fce  discusse 
briefly. 

1.  Acoustic  Waves  ir.  a  Bulk  Isotropic  y.ediun 

Jones  ana  Kobext~  have  approximately  solved  the  non¬ 
linear  equations  of  motion,  Eqs.II-5,  for  a  bulk  medium 
using  the  perturbative  approach  and  a  Green's  function 
technique.  They  take  for  U 


=  A  cos i 
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anc  t  r  i  s  coiTiGiricitior  01  two  prirr.ary  v'aves  geners,  ^6£  t .r.e 
three -phonon  interaction . 

From  the  homogeneous  solutions ,  which  are  both  shear 

and  longitudinal  plane  waves,  Jones  and  Kobett  construct  a 

— * 

Green’s  function  G  (r,  r’,<u)  for  Eq.II-5.  The  solution, 
after  taking  the  time  Fourier  transform  of  Eq.II-5,  in  the 
infinite  region  can  be  written  as 


U'(r.ca)  =  fv  G(r,r’,cii)  q(r',o>)  av 
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^where  G  is  a  tensor  of  rank  two  ana  consists  of 
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is  expresses  as 


The  solution  is  quite  complicated  and  the  original 
paper  should  ce  consulted  for  details.  However,  they  fi 
xhat  under  a  certain  resonance  condition  the  solution  U’ 
an  independently  propagating  wave.  Hence  a  third-phonon 
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The  resonance  condition  also  gives  the  limit  to  the 
ratio  of  the  primary  wave  frequencies.  In  this  limit  we 
can  choose  an  angle  between  the  primary  waves  vectors  k^ 


an: 


:c>  that  we  ??:  a  ffeneratec  wave  m  tne  a  ore' 
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Another  important  result  is  that  the  amplitude  of  U’  is 
proportional  to  the  product  of  the  primary  wave  ampli¬ 
tudes  . 
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The  line 

ar  equations  have  been 

solved  by 

’  Farnellw  for 

an  acoustic 

wave  propagating  on  the 

surface 

of  single  cry 

s  t  a  i  s  an  p  _  e  . 
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Tanski  ha'  suggestec 

several 

pcssitle  appr 

aches  for  calculating  the  magnitudes  of  surface  wave  inter¬ 
actions.  The  equations  of  motion  given  by  He. II -2  also 
describe  the  nonlinear  interaction  of  surface  waves,  pro¬ 
vided  the  solutions  satisfy  a  stress-free  boundary  at 

x  „  =  0  , 
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Several  approaches  may  be  taken  for  solving  Eqs.II-2 

and  11-11  which  describe  the  nonlinear  interaction  of 

crystal  surface  waves.  The  first  two  approaches  suggested 
.  2 

by  Tanski  both  follow  the  perturbative  metre d  usee  tr.  tr.e 
isotropic  bulk  medium  case.  The  third  approach  assumes 
that  the  full  solution  is  composed  of  three  surface  waves , 
two  primary  and  the  resultant.  The  solution  is  constructed 
such  that  the  conservation  criteria  or  resonance  condi¬ 
tions  are  satisfied.  This  gives  the  important  result  that 


l^the  acoustic  wave  interaction  for  plane  waves  in  an 
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.  Solving  the  Acoustic  Equations  for  Harmonic  l-er.eraticn 


acoustic  v;ave  of  a  given  frequency 


;,r.er.  a  sinusoidal  acoustic  wave  u:  <&  s- 
ani  sufficient  amplitude  is  introduced  into  a  nonlinear 
cr  anharmonic  solid,  the  fundamental  wave  distorts  as  it 
propagates,  so  that  second  and  higher  harmonics  of  the 
fundamental  frequency  will  be  generated.  In  this  section 
second- harm o n i c  generation  in  an  isotropic  ouik  men um  and 
crystalline  surface  waves  will  be  discussed  briefly. 

1.  Acoustic  Second -harmonic  Generation  m  an  Isotropic 
Bulk  Medium 

The  nonlinear  ecuations  for  an  isotropic  ouik  solid 
are  g  —  —  n  c  j  a  q  •  -  x  —  o  •  1  up  p  c  s  e  w  d  ^  ^ 

1  longitudinal  wave )  is  introduced  into  an  isotrctic  bulk 

medium.  The  amplitude  th  may  be  written  as 


U,  =  Aq  sin(k1x-o)-L‘t ) 


and  in  this  case  Eq.II-3  reduces  to 
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An  analysis  of  second-harmonic  generation  in  an 
Lisotropic  bulk  medium  reveals  three  important  results: 


J 


i).  A  necessary  condition  for  acoustic  second -harmonic 
generation  .s  that  no  dispersion  cf  chase  velocity  is 
carries  idle .  The  velocities  cf  the  primary  and  secondary 
waves  must  be  the  same;  that  is 


V_(  2cs.  )  =  V,  f<o,  )  . 
a  1  1  l 


11-15 


where  is  the  velocity  of  the  second  harmonic.  (2).  If 
a  secondary  plane  wave  is  generated  and  ,  then  the 

amplitude  of  the  second-harmonic  wave  is 

A„  a  A*T  x  .  11-16 

^  0 

'3  )  ,  zf  si  secondary  nlane  wave  dees  nor  need  xr.e  recess  an 
condition,  (Ec.Il-15).  "he  amplitude  will  remain  small. 

2.  Acoustic  Surface  Wave  Second-harmonic  Generation  in 
Crystals 

The  linear  equation,  Eq.II-2,  with  right  hand  side 

terms  neglected  is  solved  for  surface  waves  in  crystals  by 

g 

Verevkina  et  al  .  The  solutions  are  subject  to  stress- 
free  surface  boundary  condition  of  an  anisotropic  medium. 
The  nonlinear  equations  of  motion  have  been  solved  appro¬ 
ximately  by  L^pen^,  using  the  perturbative  approach  dis¬ 
cussed  in  Sec. II  B.l.  From  the  solution  we  know  that  the 
amplitude  of  the  second-harmonic  displacement  vector  is 
.proportional  to  the  interaction  length  and  the  square  of 
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CHAPTER 


ISOTROPIC  FLATS 


In  the  preceding  chapter,  we  have  reviewed  the  acous¬ 
tic  nonlinear  interactions  in  a  bulk  medium  and  on  a 
crystal  surface.  We  will  devote  this  chapter  to  solving 
the  nonlinear  equations  of  motion  for  an  isotropic  solid 
plate,  a  more  complicated  case,  through  a  perturbative 
approximation  method  and  Green's  function  Technique. 

The  basic  equation  cf  metier  for  the  isotropic  plate 

is  given  by  Zq.II-3-  In  practical  applications  the  dis- 

— > 

placement  vector,  U,  is  small  and  the  nonlinear  terms  m 
U  will  be  considerably  smaller.  Thus,  a  perturbative 
approach  may  be  used  to  solv  Eq.II-3*  We  set 
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III -1 


where  U  is  the  solution  with  the  nonlinear  terms  or  right 
o 

hand  side  of  Eq.II-3  set  equal  to  zero.  U'  is  considered 
a  small  correction  due  to  the  right  hand  side  of  Eq.II-3 • 
Both  UQ  and  U'  are  subject  to  boundary  constraints. 
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where  A  and  B  are  arbitarv  constants.  c  is  given  by 
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/here  k^  is  the  bulk  longitudinal  wave 
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where  k.  is  the  bulk  shear  wave  number •  The  k  are 

t  ns ,na 

the  wave  numbers  of  the  symmetric  and  antisymmetric  Lamb 
modes,  respectively.  For  an  isotropic  plate,  the  solu¬ 
tions  to  Eq.II-3  are  subject  to  two  ooundary  constraints; 
both  plate  surfaces  are  to  be  stress-free.  The  stress 
conditions  are 
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on  the  planes  z  =  "d.  In  the  linear  approximation  these 
stress-free  conditions  reduce  to 
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The  particular  modes  that  can  exist  at  a  frequency  f  ir.  a 
plate  of  given  thickness  a  are  determined  by  these  boundary 
conditions.  In  Fig.  3*  s.  graph  of  variations  in  phase 
velocity  as  a  function  of  fd  for  symmetric  and  antisymme¬ 
tric  Lamb  modes  in  vaccum  for  an  unloaded  brass  plate  is 
shown. 
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Now  the  source  tern  P-  is  comocsed  of  a  sun  of 

10 

products  of  the  two  primary  waves  in  the  three -phonon 

interaction  problem.  In  the  second-harmonic  gene rati 

case,  P.  is  comoosed  of  a  sum  of  oroducts  of  the  ini 
10 

wave  with  itself.  The  inhomogeneous  plate  equations, 
111-4,  may  be  written  in  vector  form  as 


?-*  —*■ 

a  U  '  (r,t) 

at2 


II 


L_ 


1 


A.,er,  -q  =  ?ic,  P: 
ve.::.::es  s.  r  6  ^ i ■. 


‘he  shear  and  longitudinal  wave 

and  C.  = J7k+^}/9 


x-r 


?< 


ref 


-  hi  J\  —  *  i  c-  ^ .  tr  o  —  * .  *  t:  r  w  Ji  -  c  -  c 


ransforr.  of  Eq .111-7  we 


rain 


r,=  }-C-X- 


-4.  -  — *  — * 

r  .  ^c:  7x7 xU  ’  (r  f2o)  =  4Trq'r  ) , 


„e  r curler 


.  ^  -  i.JO  d_ 


f  ^  ) 


2~r 


(r,cc;  a: 


Using  a  Green's  function  method  the  solution  to 


.11 -C 


may  be  expressed 


as 


Here  the  form  of  Green’s  function  is  dependent  upon  the 
boundary  conditions.  Two  possible  methods  for  construe 
ting  Green's  function  are  an  eigenfunction  expansion  an 
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G  for  thin  plates  using  the  image  technique.  The  main 
contribution  of  this  thesis  is  to  construct  a  Green's 
function  for  an  arbitrary  plate  by  using  the  eigenfunction 
expansion  method. 


L_ 


J 


M. 
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rc.  Green's  Function  Technique  by  E irenfunction  Expansion  ' 
Representation 

As  a  brief  introduction  to  the  eigenfunction  expansion 
representation  of  the  Green's  function,  consider  0,  a 
linear  differential  operator  with  a  complete  orthonormal 
set  of  eigenfunctions,  {y*n(r)}  ,  where 

0%(f)  -  *nyn<?>  . 

Suppose  we  are  given  the  inhomogeneous  equation 


A  "  ‘  ^  -  V  — ^ 

C  v(r)  v(r)  =  w(r), 


III-10 


where  v  is  a  vector  to  be  determined,  w  is  a  known  vector, 
and  \  is  a  constant  not  equal  to  >  for  any  n.  V.'ritting 


E  v  y 
n  n  •'n 


and 


A  — v 

0  v 


=  I 
n 


vn  Xn  yn 


equation  III-10  becomes 


z  vn  (  xn-  x  )  yn 

n 


=  w 


Taking  the  scalar  product  of  both  sides  with  y  ,  we  obtain 


L_ 


J 


J 


-  c  /  - 


(  A  -  >.  )  = 


'•  r  '  )  *  -v  ( r 


Tn is  gives 


J1  ym(?’ )  *”w(7’ )  dr’' 

5T^ X 

m 


v  = 
m 


anc 


v  = 


L. 

n 


yn(r)  J  yn(r')*w(r')  dr 


^  n~  X 


,  -  yn(r)  yn(r,)'w(r,)  dr> 


n 


\  -  k 


j  G  ( r » r  ;  *  vv ( r  ’ ;  dr  ’ 


n 


Thus  the  Green’s  function  for  the  differential  equation 
III-10  is 


?(?,?* ) 


=  z. 

xA 


n 


—  /-* ,  — ►*  ,-*■  , 

Y  (  r  j  V  (  r*  '  J 

J  n  '  f  n  ;  ~ 

*  n 


a  -  • 


It  is  important  to  note  that  a  necessary  condition  for  the 
eigenfunction  expansion  representation  of  the  Green’s 
function  is  that  the  set  of  eigenfunctions  {ya)  he  both 
orthogonal  and  complete . 


L_ 


J 


c  the  homo- 


"1 


?c r  a r.  isctrctic  glace,  the  solutions  to  tr.e  nomo- 
ger.ecjs  differencial  equations  of  motion,  Eq.III-6  with 
q-.r,--'  set  to  zero,  with  stress-free  boundary  conditions 
are  the  Lam.c  modes  of  Zc .111-2 .  To  obtain  a  Green's 
function,  Eq. II 1-11,  for  Sc. Ill one  needs  to  consider 

the  properties  of  orthonormality  and  completeness  of 

„  — ^ . 

these  Lamb  modes,  iLrc*  Both  these  properties  present 
difficulties.  The  key  problem  is  that  the  Lamb  modes  are 
r.ot  orthogonal.  However,  an  orthogonality  relation  for 
the  Lamb  modes  has  been  obtained.  A  derivation  of  this 
result  is  given  in  Section  D. 

The  seccr.i  difficult'.'  is  a  serious  one  and  vet  to  be 


:o_veo.  :cr  stress-tre- 


:cn  s  trai: 


completeness  tas 


be  assumed  and  is  not  yet  provable.  Various  authors*- 
have  in  the  past  assumed  completeness  of  the  eigenfun¬ 
ctions  (plate  modes)  with  very  good  results  which  are  in 
agreement  with  experiments.  Therefore,  completeness  will 
aloe  be  assumed  in  this  work.  A  discussion  of  complete¬ 
ness  is  given  in  Section  D  also. 


L 


J 


of  the  ' 


r. 


ni£r0n^'uinc"tiior'  G  r  t  h  c  a  o  n  ^  ^  ^  t  ^  q  y*'1  3.v'"'  ^or^'c'"  c  *  c  p.  ^  s  s 


x^cjnb  wooes 

In  this  section,  following  £  derivation  of  Fraser’s^4, 
we  will  find  orthogonality  relations  and  discuss  complete¬ 
ness  for  the  Lamb  inodes  with  various  boundary  conditions. 

i  5  ... 

maser,  using  earlier  work  by  Fama  on  elastic  vibration 
of  a  circular  cylinder,  obtained  these  orthogonality  rela¬ 
tions  by  considering  certain  adjoint  differential  equa¬ 
tions.  The  properties  of  such  equations  have  been  studied 
in  detail  by  Naimark^' .  The  orthogonalized  plate  eigen- 


■f  ijj  m 

c tic ns  which  we  wil 

1  thus 

ott 

o  '  n 

an  a 

s  c 

iutions  of 

the 

h  c  7. 

cgeneous  equat ions 

III-lu 

and 

T  — 

*"  -O 

an 

d  v;  ill  c  6  u 

sec 

to 

construct  a  Ireen’s 

function 

for 

the 

~  Y~ 

homogeneous 

ncn 

linear  equation  II-3.  First  we  must  cast  the  plate  equa¬ 
tions  in  a  suitable  eigenfunction  form. 

1.  The  Eigenvalue  Equations 

Consider  a  plate  bounded  by  surfaces  z  =t  a,  and 
assume  the  motions  are  independent  of  the  y  coordinate 
and  that  the  y  component  of  the  displacement  is  equal  to 
zero.  Let  “  and  ~  be  dimensional  variables,  we  introduce 
dimensionless  variables 


(u.v: ) 


d 


'  x  ,  z  ) 


,  and 


( U  ,  W  ) 
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where  is  the  velocity  of  dilatation  waves  in  the'  infi¬ 
nite  medium,  h  is  the  average  displacement  amplitude,  and 
U,  W  are  displacements  in  x,  z  coordinates,  respectively. 
We  may  also  expand  all  dependent  variables  in  the  form 


U (x ,  z,  t) 
W(x,  z,  t) 


V*> 

-  z  a 

n  n  W  (z) 


i(k  x-at) 
e  v  n 


III-13 


where  ca  is  the  angular  frequency  of  vibration  of  the  plate , 

and  the  wave  number  k  is  the  eigenvalue  for  the  eigenfunc- 

n  w 

tions  U„ ' z ) ,  W(z).  Now  we  may  derive  the  adjoint  differ¬ 
ential  equations  to  obtain  the  orthogonality  relation  for 
the  Lamb  modes. 

The  homogeneous  displacement  equations  of  motion,  Eq. 
II-3  with  right  hand  side  equal  to  zero,  in  components 
form  are : 
x  component: 


J 
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he  wish  to  consider  the  vibrations  of  the  plate  as 
an  eigenvalue  prctl er. .  However  the  difficulty  with  Eos. 
III-17  and  III-lS  is  that  both  k„  and  k^  appear  ir.  Ess. 
III-17,  III-19.  and  III-20.  To  get  around  this  difficulty 
instead  using  Un  and  W  as  dependent  variables,  Fraser*^ 


chose  either 

W  and  A 

,  or  U 

and  .7  . 

He  thus 

obtains 

n  n 

n 

n 

two  pairs  of 

equations 

in  which  the  e 

igenvaiue 

linearly  and 

further , 

these  equations 

are  ad  jo 

mt  xo  6  s. c r - 

other.  Take  Eqs .ill -16  ana  III-21  as  they  stand  as  the 

(W  ,  A  )  equations.  Using  Eq.III-20  to  eliminate 

from  Eq.III-19  then  we  obtain  an  expression  for  A  in 

terms  of  n  and  U  ,  which,  when  substituted  into  Ec.III- 
n  n 

17  gives  us 


L 


J 


(I-c- 


III -23 


The  (U  ,  ~  )  equations  are  given  by  Eqs.III-22 
In  order  to  express  Eqs.III-lS,  22  and  Eqs . 
in  eigenvalue  form,  we  define  the  following  two 
vectors : 


and  III-23 
III -21 ,  23 
component 

III -2^ 


HI-25 


III-26 


t  35 


Eqs.III-18,  21,  22,  and  23  may  then  be  compactly  expressed 
as 

III -2? 
III -28 

These  equatictns  are  adjoints  of  one  another  according  to 

n 

the  definition  of  Naimark  . 


A  — *>  p 

L(Rn}  =  kn  Rn 


/V  — ^  9  — 

M(Q  )  =  k  Q 
n  n  n 


L 


J 
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n, 


2.  Orthogonality  Relations  for  Different  Boundary 
Conditions 

An  equation  useful  in  considering  orthogonality 
relations  can  he  obtain  by  multiplying  Eq.III-27  by  Q 
and  Eq.III-28  by  R  . 


1 


n 


(km  -  kn>  d2  =  '-l  f«V  ' V  dz 


=  [boundary  terns  J  d 


III-29 


Here 


*  w  n  o  A  S  W  S  Q 

[boundary  terms  J  =  (1-c  )  n„  <—2  >  +  nn  ^ 


3 


au.  a 


_T]  £ —  ( — £i£ )  +  _ n  (-H) 

Un5z  K  2}  az  1  2 
c  c 


are  easily  obtained  by  partial  integration.  When  these 
boundary  terms  vanish,  a  conventional  orthogonality  re¬ 
lation  results 

III-30 

S}  "rI^L  dz  (W  o  -  ~  A  U  )  dz  =  0 

-1  m  n  -1  m  n  m  n 

c 

for  n/m 


To  obtain  an  orthogonality  relation  for  the  stress- 
free  plate,  of  particular  interest  here,  it  is  useful  to 

L.  J 


To  find  an  orthogonality  relation  in  terms  of  stress,  we 
use  Eqs.III-31  and  III-35  to  eliminate  and  r<n  from 
Eq . IT  I -30. 

L 


witr.  cc 


£  r.  o  e  s 


=  o 


;  cr 


III-3 


whenever 


111-3 


cf 


■-  '■  r*  '  —C-  ' 


ation  for  the  stress-free  plate, 

12 

Fraser  considers  four  sets  of  possible  homogeneous 


boundary  conditions  at  z  =_i 

for 

the 

plate : 

(1)  u=w=c, 

(2)  U=C  =0,  (3)  W=d  =0  and 
Z  Z  X  z 

(4) 

a 

zz 

~  ^  —  \J  • 

A  Z 

It  is  the 

last  case,  the  vanishing  cf  both  stress  at  the  surface, 
which  is  of  particular  interest  in  this  investigation. 
Although  cases  (l)-(3)  will  not  be  used  in  our  calcu¬ 
lations,  we  will  carry  through  the  orthogonality  argu¬ 
ment  for  them  in  anticipation  of  possible  future  need. 

In  cases  (i)-(3)»  it  is  easy  to  show  that 


0 


[ 


boundary  t 


-errr.sj  : 


so  that  the  conventional  crt'ncgcnality  relation  Zc_ .  1 1 1 -3  C 
hold.  The  constrains  for  the  Lamb  nodes  are  stress-free 
boundaries,  case  (4),  and  it  is  straightforward  to  show 
Zq. HI-39  holds  in  this  instance.  Thus  Ec. 111-35  is  the 
orthogonality  relation  for  the  lane  nodes. 

3-  Completeness  with  Different  Boundary  Conditions 

Casting  our  original  pair  of  ecuatior.s  Eos.  Ill-Id  and 


HI-13  into  the  adjoint  eigenvalue  form  Eqs.  HI-27  and 
III-2S  accomplishes  several  things.  A  set  of  eigenfunc¬ 
tion  solutions  is  found  and  orthogonality  relations,  Eqs. 


HI-3C  and  III-30,  for  them  are  proved.  Secondly,  system 

of  such  adjoint  equations  have  been  extensively  investi- 

13 

gated  by  Naimark  .  In  particular,  completeness  of  these 
sets  under  various  boundary  conditions  has  been  studied. 
In  case  (2)  and  (3)  the  boundary  conditions  are 

Theorem  Sec. 2.1  can  be  invoked  to  assert  that  form 


As 


L(R ) 


0 


2 

ensures  that  the  eigenvalues  are  denumerable,  and  the 
system  of  eigenfunctions  {Rn3  is  complete.  The  compiet- 
ness  theorem  clearly  holds  also  for  the  vectors  Q  . 
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^~In  case  (1)  the  boundary  conditions  are  not  regular,  and  ^ 

in  case  (4)  the  boundary  conditions  contain  the  eigenvalue 

2  16 
k  so  that  Naimark’s  theorem  is  not  relevant.  Fama  who 
n 

has  applied  Naimark's  theorem  to  thr  orthogonalized' 
solutions  for  a  vibrating  cylinder  believes  that  for  these 
cases  the  conditions  of  the  theorem  may  be  too  strong. 
Because  a  number  of  calculations  assuming  completeness 
in  case  (4)  have  yield  good  results  for  the  cylinder, 

Fama  conjectures  -that  this  set  of  solutions  is  essentially 
complete.  Based  on  close  analogy  between  the  plate  and 
xhe  cylinder,  we  too  will  assume  completeness  for  the  case 
(4)  boundary  conditions. 

4.  Expansion  of  an  Arbitrary  Function 

In  Section  III-C,  we  briefly  summarized  the  use  of 
a  complete  set  of  orthogonal  eigenfunction  in  construct¬ 
ing,  through  the  Green's  function,  a  formal  solution  to 
an  inhomogeneous  differential  equation.  Essential  to  that 
argument  is  the  expansion  of  an  arbitrary  function  in  terms 
of  the  set  of  eigenfunctions .  The  property  of  completeness 
insures  the  existence  of  the  expansion  and  orthogonality 
is  necessary  to  find  the  coefficients  in  that  expansion. 

In  the  preceding  sections  we  have  derived  from  the  Lamb 
modes  a  new  set  of  solutions  for  the  stress-free  plate. 

We  have  developed  orthogonal ity  relations  for  this  set 
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r  n 

and  have  argued  that  they  are  complete.  In  this  section 
we  will  see  how  to  expand  an  arbitrary  function  in  terms 
of  this  set. 

From  Eq.III-24,  the  new  plate  functions  are 


R  = 
n 


’  Qn 


n  q,  _ 

n  In 


Un  J  Q2n 


II 1-40 


in  which  A  and  O  are  given  in  terms  of  U  and  V!  by 
nn  n  n 

Eqs.III-19  and  III-20  respectively.  The  orthogonality 
relation  for  cases  (l)-(3)  is 


-r.i  *.fi,  <*vv  4  - 


=J'-l  <RlmQln  *  Q2nR2m>  d2  =  0  if  m  /  n  II I -41 


and  for  case  (4) 


,  ,  P  au 

r1  (w  o  -a  u  )  dz  =  rJ-  { w  c  (n  +2  — -)  - 

J  _1  Vxzn  xxm  n'  J -l**  m  n  az 


<%-2d' 


dW  9  -i 

-jf>  »„)  dz  -  c  J~-l  £  Rlrl«ln"2  ^  *  <R,.»  + 


2  SRlm 


)  Q?  }  dz  =0  if  m  ^  n 


dz  '  2n 
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E .  Formal  Solution  of  the  Inhomoe'ere: 


.e  ucuations 


n 


In  this  section,  we  will  use  the  eigenfunctions 


n 


i.na 


ln  to  solve  the  inhomogeneous  nonlinear  equations., 
Zqs.III-6,  for  an  isotropic  plate.  In  component  form,  we 
may  write  Eq.III-o  as  x  component: 


o  a  *"u 1  ,d2u'  _  a2u' 

>  o  -?  '  ^-TT"  + 


c  X 


c  2 


CX 


+  T=pr )  =  P  111-4-6 
0Z0X  a  0 


and  z  component: 


? 


o 


a 


) 
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In  Section  D,  we  manipulated  the  homogeneous  versions  of 
these  equations,  Ecs.III-i4  and  Eqs.III-15,  into  the 
homogeneous  adjoint  form  Eqs.III-27,  III -28.  In  this 
section,  after  dealing  with  the  "x  and  t  dependent  vari¬ 
ables,  we  will  carry  through  essentially  the  same  manipu¬ 
lations  on  Eqs.III-46  and  111-47  to  obtain  an  analogous 
inhomogeneous  adjoint  form. 

If  we  assume  an  exponential  dependence  in  x  ana  t 
and  take 


U  '  (x , z , t ) 


U  ’  (z) 


i (k ’ cos  G ' x-v 


t ) 


L_ 
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n 


“.'-T-'i  -  “ .  /  „  N  _i(k  ’cos&x-i  '  t) 

Equation  III -46  becomes 
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?0  CJ  (U‘  ( iio  '  )23  -  u{U' (ik-cose- )2  +  ^2  3-  (K+|) 

o  Z 

i  ^i(k'cose’x-co't) 


it."  (ik 'cosS' )2  +  'cos  6’ ) 1  e 


fxc(x,z,t) 
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Later  when  we  examine  the  source  function  P  (x.z.t)  we 

xo 

will  fir.i  ir  to  have  ar.  x,  t  cecer.ier.ce  of  the  form 


r  w  ;  e 

xo  ' 


and  the  validity  of  Eq.III-49  at  all  x  anat  will  imply 
that 

k'cos©'  =  k  ,  CD"  =  CD  ' 


Thus  one  exponentials  can  be  divided  out,  and  we  are  le: 
with  only  the  z-dependence  of  Eq. 111-^9* 

In  anlogy  to  Eqs.III-19  and  III-20,  we  define 


A  '  =  ik  ’U  ' 

+ 

a  z 

III-50 

r. '  =  ik  'Vv 

-  -  1 

0  u 

III -51 

'""in  terras  of  which,  after  some  manipulation,  Eq.  111-49 
becomes 

(  2  „  p  • 

(1-c2)  +  Mil  +  W,2U*  -k,^cos2fe'U’  - - III-52 

oZ  ^  2  pr2 

o  Z  y  C 

o  Z 

Similarly,  Eq.III-47  can  be  written  as 


i-c^  |a_L  +  ajr.  +  _  k,2cos2e*w  = 

c  az  c 


hci=2 


III-53 


Also,  in  analogy  to  Eqs.III-21  and  III-22,  we  can  obtain 
equations  for  n’  and  '  individually: 


2  ^  dP  * 

*  <4"-  fc'2c°s2e'>  n'  *  ~j(-  -sr  -ik'cos6T'o> 
a*  c  fo  ct 

111-54 


MU  +(  a>,2-k,2cos2e') 


A’  =  — U(-  t22—  -ik  'cose  'P '  ) 

f  ^2V  az  xo 
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In  order  to  express  Eqs . III-53 ,55  and  Eqs.III-52,  54 
in  eigenvalue  form  we  define 


a-c2)  k'n?)  -  *  S'"' 

L(R’)  =  c  az  c 


,2  ,  ,2 
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o  Z 


III-57 


Then  These  pairs  of  equations  may,  respectively,  be  writ" 
en  as 


L(R’)  -k ’ 2cos2e  *R  ’ 
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+  ik *  cos  0  *  P  * 


xo 


=  P- 


III-58 
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-lk’cosO’P’ 
s:  20 
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-P’ 

xo 


=  p. 


z1 


II 
III-59 


For  convenience  and  to  emphasize  the  analogy  with  Eqs . 
111-27  and  III-2S,  we  denote  the  source  vectors  on  the 
right  side  by  p ^  and  P-^  so  that  Eqs. 111-58  and  III-59 
become  finally 
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L(R  ’  )  -k’2cos2e,R'  =  P 


III-60 
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M(Q ' )  -k '  cos  0 ' Q '  =  P 


II 
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Considering  Eqs.III-2' 
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K(Q  )  -k  Q  =  0 
v  n  n  n 


III -25 


and  Ecs.III-60,  6l  together,  we  may  now  easily  obtain  a 
series  solution  to  the  inhomogeneous  equation  following 
the  procedure  sketched  in  Section  III. 5.  1/ve  have  argued 
in  Section  III.D.3  that  the  function  R  are  complete  so 


that  both 


the  given  source  function  P.  ad  the  desere: 


solution  R’  may  be  expanded  in  terms  of  them. 
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Inserting  Eas.III-62  into  III-60 
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n  n 


ana  eliminating  L(F>n)  through  Eq.III-27,  yields 
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for  the  known  function  F_  may  be 


detained  frcir.  Ec.III-bh  when  the  boundary  conditions 
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Finally,  we  have  the  solution  R'  of  the  inhomogeneous 
Eq.III-60  as 
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It  is  straightforward  to  show  in  terms  of 
the  appropriate  Green’s  functions  are 
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and  ,  finally,  the  solutions  cf  the  inhomogeneous  Zqs 
III-60  and  III -61  for  stress-free  boundaries  are 
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chapter  IV 

CALCULATIONS  AND  RESULTS 

In  the  preceding  chapter,  we  formed  a  formal  solution 
for  the  nonlinear  equations  describing  the  vibrations  of 
a  homogeneous  plate.  This  solution  takes  the  form  of  an 
infinite  series  Eqs.III-66  and  111-6? »  and  may  also  be 
expressed  in  Green's  function  form  Eqs.iii-7?  and  III-75. 
Formal,  infinite  series  solutions  are  pratically  useful 
either  when  they  can  be  reduced  to  'use  a  few  terms  for 
calcuiat ional  purposes  or  when  some  general  property  can 
be  discussed  by  the  examination  cf  a  single  term,  he  will 
consider  both  these  possibilities  for  nonlinear  plate 
solution  formed  in  Chapter  III. 

Although  an  extensive  numerical  calculation  is  now 
possible,  we  will  restrict  ourselves  here  to  a  consider¬ 
ation  of  the  symmetry  properties  of  the  solutions  arising 
in  the  three-phonon  process  and  in  harmonic  generation. 

Experimental  measurements  on  symmetry  properties  with  which 

3 

we  can  compare  our  results  have  been  taken  by  Brower  . 

Under  certain  experimental  conditions,  a  mechanism 
exists  for  isolating  a  single  term  or  limited  group  of 
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1  terms 
last  c; 
forms  f 
factor 


n 

our  infinite  ser.es  plate  solutions  formed  m  tne 
ter .  Fror.  either  the  series  cr  Ireer. ' s  fume t ton 
is  seer,  that  the  series  coefficients  certain  a 


(  k2  -  k  '  2cos2c '  )_1 


iv-: 


in  which  k’  and  cos  O'  are  determined  by  the  momentum  and 
relative  angle  of  the  two  input  pump  modes.  kn  is  the 
magnitude  of  one  of  the  allowed  modes  of  the  plate  as  dis¬ 
cussed  in  Sec.  II. A  and  illustrated  in  Fig.  3-  Thus  when 
the  combined  momentum  of  the  pump  modes  is  close  to  that 
o i  an  a _ owes  p_l ate  mcce  so  that 

2  2  a*-  2 

k  '  cos^e-  = 

n 

and  the  factor  (Ea.IV-1)  is  large,  the  nth  term  can  be 
expected  to  dominate  the  series. 

In  this  chapter,  we  will  study  the  properties  of  a 
single  term  in  the  series  solution  with  two  objectes  in 
mind.  First  to  determine  the  symmetry  characteristics  of 
this  term  associates  with  the  pump  modes.  These  results 
are  of  immediate  interest  in  the  present  investigation. 

And  secondary  to  simplify  and  partially  evaluate  the  term 
in  anticipation  of  future  numerical  evaluation. 

In  the  following  sections  five  cases  will  be  studied. 

L  J 
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These  are  the  three -phcr.cn  process  for  combinations  of  the 
two  pump  modes  which  are  antiaymnetric-antisymme trie ,  anti¬ 
symmetric  -syrnmetri  c,  symmetric-symmetric  and  harmonic 
generation  for  symmetric  and  antisymmetric  pump  modes. 
Noting  that  these  processes  are  all  generally  quite  similar 
we  will  treat  the  first  two  in  great  detail.  The  last 
three-phonon  case  followes  easily.  And  recalling  that 
harmonic  generation  can  be  regarded  as  a  special  case  of 
three-phonon  interaction  with  coilinear  pump  modes  of  the 
same  frequency,  results  for  these  processes  can  be  written 
dewn  without  much  additional  work. 

A .  Three -Phonon  Interaccion  in  an  Isccroric  Plate 

To  investigate  the  three-phonon  portion  cf  the  inter - 
action,  two  primary  waves  must  be  present.  Hence  U  des¬ 
cribed  in  Chapter  III  is  now  composed  of  two  initial  Lamb 
modes . 


Uo  = 


Ao  L1  + 


B  L0 
o  2 


V-2 


where  A  and  B  are  constants,  and  (n=l,2)  are  sciu- 
oo  n 

tions  to  the  linear  homogeneous  equations  of  motion  with 
the  stress-free  boundary  conditions. 

In  the  following  subsections,  we  will  study  three 


combinations  of  the  two  pump  modes  bases  on  their 
symmetry  characteristics . 
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1.  Three-phonon  Interaction  of  Two  A 
Modes 

For  the  three -phonon  interaction 
Lamb  modes ,  UQ  is  expressed  as 
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Disregarding  "self-interaction"  terms,  terms  with  A^  or 
BE ,  PQ  may  be  written  in  component  form,  with  all  sub¬ 
scripts  a  omitted. 
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And  for  stress-free  boundaries,  case  4,  the  orthogonality 
relation  is  Ec  .III  —^-7 


(  I  •  J  =0  for  n  /  m 
n  in 

Also  we  select  the  appropriate  Green's  functions,  Ecs. HI- 
75  and  III-76. 


?■ 

G  (  z  ,  2  '  ) 


Rn(a)1;(z') 


T^(k^-k’^cos^6’ ) 


H(z  ,z  ’  ) 


XU,~4U,) 

T'(k_-lt  ,2cos28'  ) 
n  n 


Since  Lamb  modes,  L  »  may  be  symmetric  or  antisymmetric 

with  respect  to  the  median  plane,  so  the  Green's  function, 

G,  is  composed  of  a  sum  of  symmetric  and  antisymmetric 

Green's  functions,  G  ,  G  .  Hence  Green’s  functions  for 

s  a 

a  stress-free  isotropic  plate  may  be  expressed  as 


L_ 


J 


-  60  - 


9 

wcos 


na 


n 


V 

1 


l 


ns(kns-k’2cos2e,) 


iv -5 


*na  'na 


na  T  ’  (k2 -k ' ^cos26 ' ) 


na  ’  na 
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Noting  that  the  orthogonality  relations  for  eigenfunctions 
with  nixed  symmetric  characteristics  have  vanished. 


.  I 

x  ns 


0  for  all  it 


These  results  assure  us  that  there  is  no  Green’s  function 
with  mixed  symmetric  characteristic.  Finally  we  may  write 
the  solutions  of  the  inhomogeneous  equations  as 
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To  evaluate  the  integral  on  the  right  hand  sides  of 


Eas.IV-7i  8,  first  we  need  consider  the  integral  relations 

(_ 
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of  functions  ( f ^ ,  g^)  with  symmetric  characteristics.  For 
two  symmetric  or  antisymmetric  functions,  we  have  the 
relations 
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and  for  a  symmetric  and  an  antisymmetric  functions ,  the 
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The  immediate  results  are  that  the  resultant  wave 
is  symmetric  with  respect  to  the  median  plane  and  bec¬ 
ause  Pqx  and  Pqz  are  all  proportional  to  the  product 

of  the  primary  wave  amplitudes.  The  solutions  are 

2  2  2 

summations  over  ns  and  have  (k  -k '  cos  6')  in  the  de- 

ns 

nominator,  if  we  consider  k'cos©'  close  to  an  allowed 

mode  k  ,  the  summations  may  be  represented  by  one  term 
ns> 

L_with  kns  =  k'cos©'.  Finally  we  may  find  by  comparing  the 
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>th  sides  of  res. 17-9,  10, 
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or  m  vector  form,  we  write 
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These  two  equivalents  are  tne  energy  and  mo: 
vat ions . 

The  momentum  conservation  also  limits  the  frequency 
range  of  the  generated  third  phonon,  tne  limits  are 
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I  6  is  the  angle  between  k,  and  k0,  and  (©->  V  )  lies  in 

this  limits,  then  the  pump  waves  will  yield  a  generated, 

or  third,  wave  in  the  direction  k,  +  k,  =  k^ ,  k.,  is  an 

i  <-  j  s 

allowed  mode  in  the  limits. 

2.  Three-phonon  Interaction  of  a  symmetric  and  .Ant  isymme  - 
trie  Lamb  Modes 

We  may  write  UQ  for  a  symmetric  and  an  antisymmetric 
Lamb  modes  three-phonon  interaction  as 
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The  process  to  fine  the  solutions  for  this  case  'symne 
tric-antisymmetric)  is  similar  to  the  previous  case 
(antisymmetric-antisymmetric),  we  need  only  write  down 
the  results  without  any  more  work. 


0,ei(k’cos©,x-cD,t)_ 


R__  ( Bhk  )  (constant  3 
na  ns 


z 

na 


T'  (k2  -fc'^ecs^e*  )  f  Ct“d‘ 
na  na  o  1 


L_ 


in  wr.ic 


r. ;  cc  ns  ;ar; 


T  '  f >  ^  -  >  *  **  c  o  s  c  r:  '  N|  ? 
na  ns.  c 


( cor.; 


>nM(c  C  -SC  )  1H  F 
-j.  *n a  na  na  na  na  oz 


C *  -c '  H’  ) 
r.a  na  na 


( k  '  P  .  +  i?  '  )■> 
nx  02  " 


i,  cons 


■  lu  'na  na*  n 


H  )  f  ik  '  q 

ns  ** 


.  —  /  -  i  _  ^  i  i 

'na  "na  na”na' 


From  above  two  equations,  we  find  tr.az  F. '  ar.d  i '  are 
antisymmetric  with  respect  to  z,  and  their  amplitudes  are 
proportional  to  the  product  of  the  primary  wave  amplitude 
also  by  comparing  the  propagating  terms  on  both  sides,  we 
obtain 


k '  =  k,  +  k. 


=  a  + 
1 


Hence  a  symmetric  and  an  antisymmetric  Lamb  modes  may  i: 
teract  to  generate  an  antisymmetric  third  phonon. 

3.  Three-phonon  Interaction  of  Two  Symmetric  Lamb  Modes 
For  the  three-phonon  interaction  of  two  symmetric 
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From  atove  two  equations,  we  find  that  F.  '  and  Q 1  are 
symmetric  with  respect  to  z,  their  amplitudes  are  propor¬ 
tional  to  the  product  of  the  primary  wave  amplitudes,  and 
the  conservation  criteria  also  hold.  Hence  the  three- 
phonon  interaction  for  two  symmetric  Lamb  modes  will 
generate  a  third  symmetric  mode. 

Now  we  may  conclude  that  the  three -phonon  interac¬ 
tion  for  case  (1)  (two  antisymmetric  pump  modes)  or  case 
(3)  (two  symmetric  pump  modes)  will  creat  a  symmetric 
resultant  wave,  and  for  case  (2)  (one  symmetric  and  one 
antisymmetric  pump  modes)  will  create  an  antisymmetric 
resultant  wave  subject  to  conservation  criteria  laws. 
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-  •  become  Harmonic  Generation  m  an  Isotropic  Plate 
To  investigate  the  second  harmonic  generation,  we 
only  cr.e  pump  mode,  hence  U  is  now  expressed  as 


n 


need 


U0  -  A  L 

where  A  is  a  constant,  and  L  is  the  Lamb  mode.  By  compa- 
risicn,  this  expression  is  a  special  case  of  Eq.IV-2,  hence 
the  second  harmonic  generation  can  be  treated  as  a  special 
case  of  three-phonon  interaction  with  collinear  pump  modes 
of  the  same  frequency. 


1.  Symmetric  Pump  Mode  Harmonic  Generation 

For  the  symmetric  Lamb  mode  second  harmonic  generation, 
is  expressed  as 
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where  A  is  a  constant  and 
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L  are  symmetric  solutions  for  the  linear  homogeneous 
equations  of  motion  with  stress-free  boundary  conditions 
at  z  =  *  1  of  an  isotropic  plate.  Considering  ’'self¬ 


interacting”  terms  only,  the  source  term  F  may  be  written 
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in  which  the  constants  d^-d^q  are  complicated  expressioi 
given  in  Appendix  II. 

The  process  to  find  the  solutions  for  the  harmonic 
generation  is  similar  to  (also  simpler  than)  the  three- 
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the  constants  C,-C,-,  are  complicated  expression  giver,  in 

i±/ 

Appendix  II. 

From  above  equations,  we  find  R '  ar.d  1'  are  symmetric 

with  respect  to  2,  and  their  amplitudes  are  proportional 

to  (the  square  of  the  initial  amplitude).  The  sclu- 

2  o 

tions  are  also  summations  over  ns  and  have  (fc  -k  ’  “* )  in 
the  denominator,  if  ;ve  consider  >*  close  to  an  allowed 
mode  k  ,  the  summations  may  both  be  reduced  to  one  term, 
with  k'  =  k  .  Finally  by  comparing  the  propagating  terms 
on  both  sides,  we  obtain 


k '  =  2k  „  and  a> '  =  2  ©  . 
ns 

Those  are  the  statements  of  energy  anc  momentum  conserva¬ 
tions,  furthermore  those  conditions  may  act  as  the  necess¬ 
ary  conditions  for  Lamb  mode  second  harmonic  generation. 

In  compact  form,  same  as  Eq.II-15* 
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2  2 
r  T  ^  ~  •  xf  T  »  TL'^r  T  •  u  - 


’  4+fh ^  + 

H2H 1  )d z  * 


1  / 


We  ray  find  fret,  above  equations  that  ’  and  Q, 1  are  symme- 

tric  with  respect  to  z,  and  their  amplitudes  are  propor- 
2 

tional  to  A  .  By  comparing  the  propagating  terms  on  both 
sides  gives 


k  ’  =  2  k  ,  and  »  •  =  2  s 
ns 

Hence  the  second  harmonic  generation  by  an  antisymmetric 
Lamb  mode  has  symmetric  characteristics. 

We  can  conclude  now  that  a  symmetric  or  an  anti¬ 


symmetric  pump  mode  will  create  a  symmetric  second- 
harmonic  but  no  asymmetric  harmonic. 
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C  .  Bounded  Purr  Vode  Approx imatior 

In  the  preceding  sections,  we  have  discussed  several 
cases  of  three -phcnon  interaction  and  second -harmonic 
generation  with  plane  pump  waves.  In  reality,  however, 
plane  waves  do  not  exist,  hence  the  three-phonon  interac¬ 
tion  or  second-harmonic  generation  is  generated  in  a 
"exciting  region"  by  bounded  pump  wave(s).  To  introduce 
the  bounded  pump  waves  in  the  nonlinear  equations  of 
motion,  we  have  to  multiply  a  function  g(x)  on  source 
term  P  .  g(x)  is  zero  outside  the  "exciting  region",  so 
it  can  modify  the  source  PQ  to  have  nonlinear  interaction 
only  within  a  region.  Thus  the  modified  eigenf unction 
equations  may  be  expressed  as 


A 

L 


’  (S’ 


)ei(k'cose'x-a>'t)  (  z  )g(x  )e 


i  f  (k^  +k2cos0  )x- 


(CDl+“2)tj 


IV-21 


M(Q '  )e 


i (k  ’  cos©  ' x-  ^t ) 


'=P^( z )g(x )ex 


^l+^2cos6 


x-  (*>-,  + 

X 


IV  -22 
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where 


?  '  ( z 


o: 


+  i k 1  cos  6 ’? 


s 


ox 


s 


P  ’  (  2  ) 
m 


»oC* 


>  '  QQC  0  f  p 

oZ  “  "  02 


C 

-P 


OX 


L ’ (R  1  )  = 


L(R*  ) -k  ’  ^cos 2  9  ’R  ' 


X\  —4 

V  '  /  o  * 


tl(Q  ’  )  -k  ’  ^cos2© ’Q ' 


For  second-harmonic  generation  case,  we  set  0'=O,  0=0, 
and  k-i^kg  in  .qs.IV-21  and  IV-22. 

Taking  the  Fourier  transform  of  Eqs.IV-21,  one  obtains 


r  l.(R>)«i(l,’eo*®">k)x*'i»’t  «X  «£  p:(2)g(x) 


ei(k1+k2cosd+k )x 


■i(^1+a>2)t 


dx 


or 
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1 


/"\ 


(?. 


- 


(k’cos6'+k)  e 


■i- 


( z )g(k1^k2co 


k  )e 


IV-23 

where  g(k^ ^k^ cos 0+k )  is  the  Fourier  transform  cf  g(x), 
and  fc(k'cosS'+k)  is  zero  when  k’cos6'/-k.  For  k/-k'ccs6’, 
we  have  /g,(k^+k2cos6+k)  =  0,  because  fc(k  ’cos6'+k)=C.  For 
k=k * ( -cos0 ' ) ,  Ea. IV-23  may  be  rewritten  as 

L •(*?*)  6  (0)e"1^''t=  pj*(z  )g(k1+k2cos6-K' cos0'  )e_i  (~1+S2  ^ 


where 

g^k^+k  cos6-k  ' cos9 '  )  =  £*  e^^l+ic2cos 


Sc'  )x  /  - 

g ( x  ;gx 

IV -24 


The  integrand  of  Ea.IV-34  oscillates  as  one  integra¬ 
tes  over  x.  The  frequency  of  oscillation  is  determined 
by  the  coefficient  of  x  in  the  exponential  term,  so  the 
frequency  is  given  by  (k^+k^cosG-k ' cos© ' ) .  The  integrand 
will  oscillate  unless 

k.,  +k2cos©  =  k’cose’ 
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on  is  called  a  resonance  cc 


rewrite  Eq . I V - 2i  as 


on. 


ns 


gr^kI+k2Cos6-“ 'cose' )=g(0)=J^ 


g ( x ) ax  =J°;  g(x)dx  . 


A/ 

Hence  gCk^+l^cosS-k  '  cOse  '  )  = 


0  k=-k  * cos6 '  IV -25 
'g/O)  k^+k2COs6=k  ' ccs@  ’ 


=  -k 

If  we  assume  g(x)  a  step  function,  then  g(  0 )=Kx ,  where  K 
is  a  constant.  We  may  then  also  rewrite  Eq.IV-25  as 


s(k. +k-cos9-k ' cos& ' )  =  Kx  t  (k. ~k~ccsS-k )  IV -26 

_  ^  j.  ^ 

Substituting  Eq.IV-26  into  Eq.IV-23  and  Taking  an  inverse 
Fourier  transform,  we  may  obtain 

vr^  !•(?•)  6  (k  '  cos©  '  +k  )  e"1^k’x+m't^dk  = 

y50  P'(z)Kx  (k,  +k,cos6+k )  e-^  ^x+  ('xi  +^2  '  0  -  ik 

-CO  f  1  £ 

orI’('R’')ei(l''cose'-“‘'t>^(z)Kxei  C  ( Vk2ctjs6)x-<ah°2,'t) 

l 

IV-2? 

Following  tne  same  procedure,  we  may  transform  Eq.IV-22 
to 
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(  z  )Kxe 


i  i 


(k-,  -'■k.-ccsG  )x- (- 

J.  c 


, - 
x  x 


r/.-2£ 

By  solving  Zqs. 17-27  and  IV-2S  through  a  Green's 
function  technique,  we  find  the  resultant  wave  is  propor¬ 
tional  ~o  the  interacting  volume  for  three-phonon  inter¬ 
action  or  propagating  distance  for  the  second-harmonic 
generation,  respectively. 

how,  we  know  that  symmetric  or  antisymmetric  Lamb  mode 
will  generate  only  symmetric  second-harmonic  wave  with 
amplitude  of  the  resultant  wave  being  proportional  to  A  , 
the  square  of  the  initial  wave  amplitude.  We  also  find 
that  two  symmetric  or  two  antisymmetric  Lamb  modes  inter¬ 
act  will  generate  a  symmetric  third  phonon,  and  one  symme¬ 
tric  and  one  antisymmetric  Lamb  mode  interact  will  give  an 
antisymmetric  third  phonon.  The  amplitude  of  the  third 
phor.or.  is  proportional  to  the  product  of  the  two  primary 
wave  amplitudes.  Ir.  reality,  an  "exciting  region" 
caused  by  bounded  pump  wave  will  give  a  volume  or  distance 
dependent  resultant  wave  amplitude. 


L_ 


^r.  ir.e  creceuur.:!  ?r. ^ r  *  f  r  ,  *6  have  € Va-Ua*6c  a  s ir^i6 
term  cf  cur  general  series  rcluci.cn  for  the  nonlinear 
vi oration?  of  the  elate  .  This  /.as  ccr.e  with  s re c is_  empr - 
as  is  on  symsne try  rrcpsrcies  and  all  possible  c c mb i rat i o ns 
cf  s yrc:ce trie- an t i s ymme trie  pump  nodes  were  considered  for 
both  the  three-phonon  interaction  and  harmonic  gene ratio, 
Cr.e  if  the  reasons  for  doing  this,  as  exrlair.ee  in  Sec* 

IV,  A,  was  to  facilitate  c  c  it  r  arisen  with  the  expe  r  inert  a- 
results  cf  Brower" . 

A .  Three -phonon  Interaction 

Two  important  results  have  been  obtained  for  the  three 
phonon  interaction  in  a  bulk  medium  or  a  single  crystal. 


There  are  the 

necessity  of 

the  conservation 

1 1 

eria 

on 6  lac  c  *na t 

nhe  generated 

wave  is  proport 

ional 

C  W' 

product  of  the 

primary  wave 

amplitudes,  i/ve 

have 

feu: 

these  results, 

Eq.II-8  and 

I I - 9 .  to  also  be 

true 

for 

Lamb  mode  three-phonon  interaction  in  the  plate. 

Brower  has  demonstrated  that  two  input  Lamb  modes 
generate  a  third  phonon  in  a  brass  plate  ar.d  that  the 


L 


harmonic  generation:', 
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5 .  Harmonic  Generation 

Studies  of  second -harmonic  generation  in  bulk  media 
and  single  crystals  have  yielded  two  important  results. 

One  is  the  necessary  condition 

V2(2cu)  =  V1(o)), 

that  is,  the  second-harmonic  must  have  the  same  velocity 
as  the  initial  pump  wave.  The  other  is  that  the  amplitude 
of  the  generated  harmonic  varies  as 

«<  A?  x 
£  1 

where  A-j  is  the  initial  pump  amplitude  and  x  is  the  dis¬ 
tance  the  primary  wave  has  traveled  in  the  medium.  The 
first  result  is  true  for  our  plate  solutions  of  Lamb  mode 
second-harmonic  generation.  The  second  was  obtained  in 
Sec.  IV. C,  where  we  modify  the  basic  plane  wave  approach 
to  take  account  of  the  bounded  volume  of  interaction. 

c 

In  Viktorov's-'  analysis,  it  is  stated  that  only  for 
the  zero  order  symmetric  pump  mode  (S  )  at  thin  plate 
frequency-thickness ,  fd  1,  is  it  possible  to  have 
second-harmonic  generation.  This  is  because  for  small 
fd's  the  velocity  in  SQ  mode  is  not  dispersive.  Later, 
however,  Brower  pointed  cut  that  harmonic  generation  in 
(_Lamb  modes  is  indeed  possible  for  particular  sets  of 


J 


_  P  c.  _ 
w 


dispersion  curves.  He  demonstrated  experimentally  that 

for  A-,  mode  of  a  brass  plate  with  thickness  b=G.S5mm  at 

f -  ?y.Hs  and  f=5i'*Kz  harmonic  generation  occurs  to  produce 

the  A0  mode.  Thus  the  S  mode  at  thin  elate  fd's  is  not 
2  o 

the  only  mode  to  permit  harmonic  generation. 

In  Sec.IV-B,  we  found  that  either  a  symmetric  or  anti¬ 
symmetric  Lamb  mode  will  generate  only  a  symmetric  second- 
harmonic,  a  result  inconsistent  with  Brower's  observations 
of  (A^»  k^)  for  the  isotropic  brass  plate.  But  by  examin¬ 
ing  the  dispersion  curves  for  the  brass  plate  carefully, 
we  find  that  in  the  region  of  interest  the  curve  of  the 
antisymmetric  and  symmetric  modes  are  very  close. 

Therefore  we  conjecture  that  what  Brower  observed  was 
actually  a  symmetric  rather  than  an  antisymmetric  second- 
harmonic  . 


L 


J 


-  86  - 

r C  .  Cor.clus  ions 

The  amplitude  of  mechanical  vibrations  in  solids  is 
r.ct  infinitesimal,  and,  therfore,  nonlinear  equations  of 
motion  are  required  to  describe  the  disturbance.  The 
principal  of  superposition  no  longer  holds  and  hence  non¬ 
linear  interacxions  may  occur.  Two  possible  such  inter¬ 
actions  in  an  isotropic  plate  are  the  three-phonon  inter¬ 
action  and  second-harmonic  generation. 

A  new  approach  is  presented  to  solve  the  nonlinear 
equations  of  motion  for  the  plate.  This  approach  invo¬ 
lves  the  combination  of  the  perturbative  and  the  Green’s 
function  techniques.  The  Green's  function  is  constructed 
by  the  eigenfunction  expansion  method  after  a  suitable 
orthogonalization  relation  for  the  Lamb  modes  has  been 
derived . 

It  is  found  that  a  second-harmonic  acoustic  wave 
generated  by  either  a  symmetric  or  antisymmetric  initial 
Lamb  mode  is  symmetric  with  respect  to  the  median  plane. 
The  amplitude  of  the  generated  wave  is  proportional  to 
the  square  of  the  pump  wave  amplitude  and  is  given  by 
a  summation  over  plate  eigenfunctions  of  terms  which 
have  denominators  of  the  form  (k  -k '  ) .  In  approxima- 

I  iO 

tion,  the  summation  may  be  well  represented  by  one  term 

when  k  =  k  '  . 
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The  three-phonon  interaction  of  two  noncoilmear 
symmetric  or  antisymmetric  Lamb  modes  which  satisfies  the 
conservation  criteria  will  produce  a  third  mode  with 
symmetric  characteristic.  Similarly,  two  primary  modes 
one  symmetric  the  other  antisymmetric ,  will  produce  an 
antisymmetric  third  wave.  The  amplitude  of  the  third  wave 
is  proportional  to  the  product  of  the  pump  mode  amplitudes 
and  is  expressed  by  a  summation  over  plate  eigenfunctions. 
When  the  total  momentum  of  the  two  inputs  waves  is  close 
to  that  of  a  plate  mode,  the  summation  may  again  be 
approximated  by  a  single  term. 

The  approach  presented  in  this  thesis  is  appicabie  not 
only  to  the  plate  with  stress-free  surfaces  but  to  certain 
other  plate  boundary  conditions  as  well.  By  finding  the 
appropriate  orxhogonali7.ation  relations  for  the  plate 
modes,  a  Green's  function  satisfying  different  boundary 
constraints  can  be  constructed.  The  third  phonon  of 
the  three-phonon  interaction  may  also  have  the  following 
conditions . 


We  will  leave  those  as  futhur  investigation. 
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APPENDIX  I 

THE  CONSTANTS  e^eg  AND  g-^gg  IN  Sec.  IV.  A.  1. 

VVa-l(d2Vdl6+d20)*q2(d3q2+dl?*d19)*qlq2(d5ql'''d6q2*d18> 
e2=d15^1(d7*d12+d13ql>tq2(dS"dll+dliq2)+qlq2<d^d9ql  + 


d10q2  ^ 


e3  =  -C2(dl  +  d2ql'fd3S2+d6S2qi+dl6ql'fd20ql)-C2S2(d5ql"d17"d19 
^'-q2<d14s2td15>-C2qi(Vdi0S24d12+d13?i)'CiS2(d'qi*d5 


td9ql+dli> 

2 

•  r  *  i  «=  * 


‘■15'-dlCi'Cisl(d2"d5q2)-Ciq2(d3q2'fd17"d19)'Clsl(d6q2tdl6 


+d18+d20) 

p  2  2 

e  -s-C.'s.  (d.  a„+d„+d„  „aT+d-  -  )-q0C,  (aa+dGs1+d1 ,  )-C,  (a,  ,s,  + 
o  1  x  ‘  2  7  x  u  -  c  x£  2  1  ~  y  1  -  -  ~  x  ^  - 

d14qI+d15> 


7-eic-(<i1+<i2s^<i3s|)*cie2s2(d5*J*<J17-d:l9)+c^c1s1((i6s‘* 


-  2, 


dl6+d20)+dl8ClC2SlS2 
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e5=C4CiC2SlS2+CiC2Sl(a?^d10S2+a12^+ClC2£2(dS+d9SiTa12 )  + 

2 

ClC2(d13Sl+di4s2  +  dl^S2+d15'> 
in  which 


=  -  ia,  b.  b„k,  k0cos6(k,  +k,,cos0 )  d~=a„b,b_ik0cosG 

i.  11212  12  2  i  1  l  2 


/i=ia^b^b2(k^+k2cos0) 

d3=ik1 a2b^b2 

5*a2blb4 

d6=a2b2b3 

„=-k, k_cos6a.  b,  b- 

7  12  4  1  ^ 

ac=-k, kpeosOa^bjb^ 

9=(-^}b2b3 

d10" ' a3 / c ^bib4 

ll=4a4/2)b2b3k^ 

d12=(-a4/2)b1b4k2c°si;e 

13= ( a^/2 )b3b4ik2cos9 

dl4s(a*/2)VVkl 

.13=(a3/2  )b3bi|(-ik1k2cose)  (k 

, +k2cos0) 

i^^=a3b2b3  ( -k^k2cos0 ) 

d-L3=-a3b1bi+kik2cosS 

l^g=  ia3b3b2+(k-^+k2cos6 ) 

d19="a6blV? 

l20’-a6b2V2cose 
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in.  which 
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i  r.  wr.icr. 

-  l-£lc3°it- 

f „ =a„b_b 

4  4  j  ~ 

f„  =  a„b„b,  (-k2cos 2e) 

3  ^  jj  t  2 

f ,  ,  =  -a-b^b,.k.  k„cos6 

15  =  'abc3D2iCi 

f4=ia..,b,  b.k0ccsS 

2 

f,_,=  -ia2b?b2k^k2cosG 

V-^aViVi00*2* 

fq=  (a^/2  )b1bi+ik2cos9 

4  =  ( a  /2  vc_b„  ’  k 

-10  lV-L2“3-li 

fll=(‘a3/'2)bl'b4iklk2cos9 

f12=(a3/2)b2b3(-ikik 

fi ^=(a^/2)b  b£ 

i'l4=(-a./2)b:b2k2 

f15=  ("a4/2)bib2r:2cOS  6 

fi6='a5b:b2k:r:7CO££ 

fl7=a5blVkl 

fl8=a5b2b3ik2cosS 

f19=a6b2b3ik2COS& 

f20=a6blb2iki 
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APPEND IX  Ii 


THE  CONSTANTS  b, -b, -  AND  C, -C 
1  lu  1  i' 


>6C  .  ±V  .  t.  1 


=  (a2+a^  )k^q2-a^k^q^+  ( a2+a,  )kq 


52=-2a1C  ’k5+(a2+a- )C  ’kJc2+a2C  ,k3s2+a,Ck-?cs+a6C  ’k2(qs2+ 


2  2 


q ^  )- ( 2ac  +  a2 )Cka^s-a2C 'kq^s 


b^  =  a^C  ’  2k3-a2C  ’  2k3s  2-a^CC  'k^as-a^C  '  ‘H-c^qs^-i-a^C  ’  2kqs 


'^k^cs^-ra.c  ,^kns'/  + 


,  2 

a-xc  s, 


Dr~  (  —  a-5"*’2  /  ^  a  i  )  k  3  g  -  ^  a  a  j  k  g  +  ( a  ^ /2  —  a  ■ ,  /  2  )  kc 
v  3  •*  j  j  -*■ 


b5=(a3-  3/2  a4)Ck3q2-(3/2)a4C ’qs+(-a3+a4)ikq2s2+ 


ta4Ckqil’-ta3C  ’kq3s 


b,  =  -a^C  ,2k3£2-ta_C2k3G2+(3/2)a„CC  'k3as+ia-CC  'kas3-ta,,C2kq2 

O  j  J)  ~  J)  ~  ^ 


b^  =  a1q3+(a2+ia3+ia4+a3)ki<'q-(a2+3/2  a3+3/2  a4+a^+a^  )k2q- 
b8=  ~  ^  a2+^a3+^a4+a5  )c  ’k\+  ( a2+a3+a4+a3  )k2q2sC+  ( ia3+£a4 ) 


8  tt2  ?a3  ?a4  a5'^  htv< 
k2s2qC ’+aAk2q3C ’ -a,q**sC 
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•V  '* 

-  cb  - 

t  =  - ( a  +*a_  jk^aC- (*a.,  +  a,  )k^sC  '  +  ( ja-  +  a^+a,  )k2a  2sC  '  + 
(a7+a?+ia,.  )Ck2c  3+ark2qs2C-a,c^s2C 

£  J  -r  O  X 

b,  c=(a9+|a;3)kUGCC  '  +  (ta4+a^)k4sC  ,2-(3/2  a4+a5+a6  )k/qs 2CC  '  - 
(a9+a^)k2a2sC2-fa^k2C ' 2+a1q2s^C2 

in  which 


a4=4a+e+2d+2b+c 


a9=2a+c 


=2a+2b 


a,.  =  2a+2a 


a,  =  tii-Q+e  a4  =  c*e 

in  which 


a=  u  +  *A 


b=K+  q/3+?A+B 


c=K-2/3  q  +  B 


d=iA+B 


e=B+2C 
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